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showed  the  need  for  more  exhaustive  testing  of  'possililc  ’'configurations  of  flap 
deflection,  slot  geometry  and  airfoil  angle  of  attack.  For  a  typical  4-element 
airfoil,  the  number  of  possible  configurations  can  easily  be  in  the  billions.  A 
new  technique,  based  on  evolution  strategy,  has  been  developed  for  the  problem  of 
optimizing  a  multi-element  airfoil  with  respect  to  its  envelope  of  raximumjJL/C  ^ 
versus  )  This  technique  was  applied  to  an  airfoil  having  a  slotted  leading  — — 
edge  flap  and  a  double-slotted  trailing  edge  flap  using  an  existing  computer 
program  for  2D  subsonic  mult i -element  airfoils  to  provide  values  of  aerodynamic 
coefficients.  Although  limitations  of  the  program  precluded  the  full  determina¬ 
tion  of  Lije  envelope  o!  (WC/-.*vs  ,  the  problem  of  maximizing  CpJ for  a 

minimum.  C.  constraint  was  solved.  The  result  of  testing  approximately  400  con- 
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over  the  initial  geotnetry.fC^he  conclusions  reached  from  this  work  are  1) 
evolution  strategy  applied  to  the  optimization  of  fnulti-eleinent  airfoils  can 
yield  substantial  improvements  in  aerodynamic  performance,  2)  the  number  of 
configurations  required  to  find  the  optimum  can  be  reduced  from  the  total 
number  possible  to  a  small  fraction  of  the  total  with  the  same  final  result, 

3)  the  inherent  simplicity  and  speed  of  the  technique  developed  lends  itself 
well  to  further  application  in  wind  tunnels,  4)  evolution  theory  appears  to 
be  the  best  choice  for  techniques  aimed  at  the  optimization  of  systems  defined 
by  a  large  number  of  degrees  of  freedom. 
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1  INTRODUCTION 

1.1  Multi-Element  Airfoil  Optimization 

It  is  no  secret  today  that  competition  among  civilian  and  •' 
military  aircraft  manufacturers  in  the  international  marketplace  is  1 
tough.  This  competition  has  demanded  improved  flight  efficiency,  and 
as  a  result  a  strong  emphasis  has  been  placed  on  design  optimization 
with  respect  to  take-off,  cruise,  and  landing  aerodynamic  character¬ 
istics  ►  Design  considerations  for  take-off  and  landing  aerodynamics 
are  primarily  directed  towards  the  reduction  of  wing  area,  increase  in 
maximum  take-off  weight,  improved  fuel  efficiency,  reduced  field  ■ 
length  requirements  and  lowered  ground  noise  levels.  The  relative 
importance  of  ’the  nondimensional  coefficients  of  lift  and  drag,  Ci 
and  Ca  respectively,  for  take-off  and  landing  are  shown  in  figure  (l)i 

Over  the  entire  regime  of  flight,  the  performance  at  take-off, 
when  aircraft  weigh*  and  engine  noise  are  at  a  maximum  for  the  flight, 
is  the  most  critical  factor  in  the  design  of  lifting  surfaces.  To 
achieve  the  maximum  rate  of  angle  of  climb,  thus  reducing  runway 
reouirements  and  noise  footprint,  maximum  Ci/O^  is  desired  for  a 
given  Ci.  The  envelope  of  values  of  maximum  Ci/Cfl  versus  Ci  for  the 
Boeing  727  wing,  obtained  from  flight  measurements,  is  given  in  figure 
(2),  and  is  more  or  less  similar  to  the  envelope  obtained  for  all 
other  airfoils2 . 

The  complexity  of  multi-element  airfoils,  as  seen  in  figure  (2) 
and  in  more  detail  in  figure  (3A),  has  been  dictated  by  the  need  for 
h  smooth,  and  xov  high-speed  transport  aircraft  shock-free  airflow 
at  cruise  conditions,  combined  with  the  high-lift  configurations  re¬ 
quired  at  take-off  and  landing2 »3.  Although  significant  effort  has 
been  devoted  to  the  study  of  advanced  high-lift  devices  such  as 
those  employing  blown  flaps,  jet  flans ,  boundary  layer  control  or 
augmentor  wings,  it  is  reasonable  to  expect  the  continued  use  of 
mechanical  devices  for  the  next  generation  of  civilian  and  military 
air c  r af  t.  Indeed,  as  shown  rn  figure  (4),  values  of  Ci  for  large, 
high-subsonic  aircraft  nearly  doubled  between  1950  and  1970  as  a 
result  of  improved  design  of  mechanical  flap  s/stems^.  Once  the 
choice  of  high-lift  system  and  associated  flap  and  airfoil  geome- 
tries  has  been  made,  the  relative  positioning  of  the  elements  must 
be  determined  such  that  the  obtainable  values  of  and  C1/C3  are 
sufficient  to  meet  performance  demands  (or  preflight  promises!). 

The  ideal  set  of  configurations  would  be  those  that  would  result  in 
the  largest  Ci/Cg-vs-Ci  envelope.  The  number  of  all  possible 
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configurations,  however,  is  very  large  even  for  the  typical  cas^3f;.'iHfa 


an  airfoil  having  a  double-slotted  trailing  edge  flap  and  -a*  slottld;  '-S.liM 

'"'-V' '  s  *  s  iVv?fJw 

leading  edge  flap,  as  depicted  in  figure  ( 3B ) . 

1..2  Empirical  Methods  -  .  •  U', 

The  most  common  means  of  determining  the  best  airfoil-flap  "{V. 
configurations  involves  many  hours  of  wind  tunnel  testing  for  a  rep-  1 
resentative  set  of  variations  of  slot  gap  width  and  overlap,,  flap 
deflection  and  angle  of  attack.  Optimum  configurations  are  then 
made  by  adding  to  the  main  airfoil  the  leading  and  trailing  edge  ’■]]. 
devices,  whose  individual  slot  geometries  have  been  optimized  with 
respect  to  only  one  adjacent  element.  An  example  of  this  empirical 
method,  reproduced  from  the  work  of  ljungstroffl5,  appears  in  figure 
(5).  Although  thin  empirical  method  can  greatly  improve  tl:e  aero¬ 
dynamic  performance  of  a  multi-element  airfoil,  there  are  two  signif¬ 
icant  disadvantages  that  must  be  considered.  First,  the  accuracy  of 
an  empirical  method  is  dependent  on  a  large  number  of  test  config-- 
urations;  however,  for  oven' the  most  well-funded  design  relatively 
few  different  eonf igurations  are  actually  tested.  For  instance,  for 
a  4-element  airfoil,  500  test  cases  is  only  5*10-6  <fa  0f  the  number 
of  possible  cases,  if  each  of  the  10  degrees  of  freedom  can  be 
varied  across  10  increments.  Increasing  the  number  of  test  cases 
comoc  at  the  expense  of  time  and  direct  operatin'?  cost  of  wind 
tunnels.  If  wind  tunnel  time  and  cost  trends  increase  at  the  rate  de¬ 
picted  in  figure  (6),  however,  it  is  doubtful  that  empirical  methods 
will  produce  better  results  than  at  the  present  time,  when  a  great 
emphasis  is  being  placed  on  improving  high-lift  device  performance6. 

The  second  disadvantage  of  empirical  methods  is  that  by  combining 
elements  optimized  with  respect  to  only  adjacent  elements,  the 
important  viscous  interaction  between  all  elements  is  unaccounted 
for.  This  interaction,  however,  is  strong  and  tests  have  shown  that 
the  combination  of  individually-optimized  components  rarely  results 
in  an  optimum  overall  eonf igurntion^ , 

1.3  Self-Optimizing  Airfoils 

It  would  be  desireable  then  to  have  a  self-optimizing  scheme, 
where  the  system  seeks  its  own  best  configuration  without  the  need  forr 
testing  every  possibili t”.  Far  from  e''ec”loJ"’ ve ,  such  schemas  hnv<» 
been  used  with  a  good  deal  of  success  in  several  different  problems 
where  systems  are  described  bv  several  degrees  of  freedom  and 
relative  'factors  of  fitness'  result  for  each  variation  of  parame¬ 
ters.  Ons  such  application  described  by  Levinsky^  attempted  to 


maximize  Cx  for  a  fixed  maximum  Cg  and  minimize  Cg  for  a  fixed 
minimum  Ci.Tor  a  flexible  2-dimensional  airfoil  mounted  in  a 


transonic  intermittant  wind  tunnel.  Results  of  this  work,  and 
results  of  further  investigation  of  a  3-dimensional  flexible  airfoil® 
in  a  transonic  continuous  wind  tunnel  are  given  in  figures  (7)  and 
(8).  In  both  cases,  the  hydraulically-actuated  leading  and  trailing 
edges  were  modified  automatically  by  on-line  computers  programmed 
with  a  gradient-strategy  optimisation  scheme.  This  is  described 
in  section  2.  The  only  human  input  was  the  factor  of  fitness  to  be 
maximized  or  minimized  (Ci,Cd}  or  volume)  and  a  constraint.  Gradient 
strategy  has  also  been  used  in  the  work  of  Hicks9,  where  2-dimen- 
sional  transonic  airfoils  were  modified  to  improve  Ci,  Cd  or 
maximize  volume.  Some  of  the  more  significant  results  of  this  work, 
performed  numerically  as  opposed  to  the  above  mentioned  wind  tunnel 


work,  are  shown  in  figure  (9).  Although  these  examples  have  resulted  1 

I 

in  substantial  improvements  in  airfoil  performance,  the  application  | 

of  their  optimization  methodology  was  found  to  be  more  cumbersome  ! 

l 

than  the  method  described  in  this  work.  Existing  optimization  | 

methods  are  compared  below  to  a  now  scheme  based  on  'evolution  1  j 

strategy'.  This  scheme  has  been  developed  to  determine  the  config-  i 

urations  of  a  4-element  airfoil  giving  maximum  Cx/Oa  for  each  value  1 


of  Cl,  where  only  the  geometry  of  each  element  is  known  initially.  f 
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2  THEORY 

2.1  Simple  Optimization  Problem-Existing  Solution  Methods 

Consider  the  optimization  problem  faced  by  an  experimentalist 
standing  above  a  room  in  which  there  is  a  single  geographical  peak 
whose  location  with  respect  to  a  2-dimensional  coordinate  system  is 
to  be  determined.  In  this  problem,  the  ’system’,  this  being  the  hu¬ 
man,  has  two  degrees  of  freedom,  his  location  with  respect  to  the 
x-axis  and  the  y-axis.  At  each  x,y  coordinate  position,  the  value 
is  determined  of  the  geographical  ’fitness  factor*,  or  the  distance 
between  the  x-y  plane  and  the  surface  below.  The  most  rudimentary 
method  to  find  the  peak  would  be  to  make  a  number  of  random  soundings 
in  order  to  obtain  a  rough  impression  of  the  surface.  In  order  to 
improve  accuracy,  the  grid  could  be  divided  into  a  fine  mesh  and  the 
height  is  then  recorded  at  each  mesh  point.  Although  this  scheme 
guarantees  that  all  possible  x,y  cobinations  have  been  checked, 
the  number  of  possible  combinations  increases  with  the  souare  of  the 
number  of  grid  divisions  and  with  a  cower  enual  to  the  number  of 
degrees  of  freedom  for  higher  order  systems.  Figure  (10)  depicts 
these  2  basic  schemes  along  with  two  gradient-based  strategies  that 
wore  developed  to  converge  to  the  peak  with  a  reduced  amount  of 
effort10. 

In  the  first  of  these  two  gradient-based  strategies,  known 
a3  the  Gauss-Seidel  Strategy,  the  experimentalist,  trying  to  find  the 
peak  with  the  least  amount  of  effort,  simply  proceeds  in  the  direc¬ 
tion  of  the  steepest  positive  gradient,  determined  at  the^  starting 
point,  until  the  gradient  becomes  zero  or  negative.  From  this 
point,  he  turns  in  the  direction  of  a  new,  locally-steenest  gradient 
and  continues  in  this  now  direction.  This  procedure  is  continued 
until  the  peak  is  found,  that  is,  a  position  has  been  reached  on  the 
x-y  plane  where  no  posit ivrt  gradients  exist  below.  The  second  of 
these  technioues,  known  as  the  ’general  gradient  strategy',  is  sim¬ 
ilar  to  the  first  except  that  at  each  stop  the  magnitude  of  gradients 
in  all  directions  is  re-evaluated.  Although  the  second,  scheme  may 
result  in  the  shorter  total  distance  between  the  starting  point  and 
the  peak  with  resnect  to  the  first  scheme,  it  reouires  more  work 
per  step.  Clearly,  the  use  of  either  of  the  gradient  strategies 
represents*,  a  substantial  improvement  over  the  inspection  of  every 
possible  combination  of  the  variable  parameters.  It  is  noted  that 
in  order  to  reduce  the  amount  of  effort  involved  in  scheme  B  of 
figure  (10),  one  is  t emoted  to  iocrers"  the  grid  spacing.  The 


danger  of  this,  of  course,  is  that  there  is  a  strong  ’pdssibllii^'^^^' 
that  the  peak  will  he  completely  missed?  this  is  precisely  the 
point  mentioned  previously  related  to  optimizing  airfoils  by  ' 

empirical  methods.  The  danger  of  using  gradient  strategies  in 
more  complicated  systems  is  that  only  a  relative  optimum  may  be 
reached,  as  the  smaller  of  two  peaks  for  the  simple  problem  just  . 
described. 

2.2  Simole  Optimization  Problem-Evolution  Stategy  Solution 

A  fifth,  and  relatively  new  optimization  scheme  that  may 
further  reduce  the  effort  reouired  to  improve  a  system  described  by 
a  number  of  degrees  of  freedom  ic  based  on  ’evolution  strategy*. 
•Relatively*  new  is  stressed,  because  the  basis  of  this  strategy  is 
derived  from  the  theories  of  natural  selection  first  postulated  by  the 
19th  century  biologist,  Charles  Darwin.  As  in  biological  systems, 
which  in  their  simplest  classifications  can  be  grouped  etccording  to 
a  number  of  parameters,  such  as  size,  number  of  appendages,  means 
of  reproduction  or  mobility,  an  engineering  system  survives  if  the 
comUnation  of  each  of  its  defining  parameters  results  in  a  ’beast* 
that  is  superior  to  all  othor  combinations.  Only  recently  has  this 
concent  been  applied  to  realistic  physical  problems,  the  majority 
of  the  work  to  date  being  attributed  to  Professor  Ingo  Rechenberg, 
of  the  Technischcn  Universitat  Berlin,  In  ord  r  to  describe  thin 
scheme,  we  consider  again  the  task  of  the  experimentalist.  Using 
evolution  strategy,  rs  depicted  in  figure  (11A),  a  number  of  test 
points  are  selected  randomly  over  a  region  near  the  starting  point. 

The  height  is  recorded  at.  each  point,  then  the  o.  pcrimentaliot 
moves  to  the  highest  of  these  points,  /ro’i  there,  the  same  number 
of  new  points  are  randomly  chosen,  and  the  procedure  continues  as 
before. 

The  distance  from  one  'hi^h  point*  to  each  new  point  in  a 
successive  set  is  controlled  bv  the  previous  distance,  that  is  if  a 
large  step  resulted  in  a  larger  height  than  a  smaller  step,  all  new 
points  would  be  determined  using  steps  approximately  eoual  to  the 
large  step.  This  important  factor  of  evolution  strategy  allows 
large  steps  to  be  taken  early  in  the  optimization  process,  and 
decreasing  step  size  as  the  optimum  position  is  approached.  The 
inclusion  of  large  mutation  steps  also  prevents  the  convergence  to 
local  maxima,  a  problem  associated  with  gradient  strategies.  In  the 
case  of  a  local  maxima  or  small  barrier  in  the  optimization  process, 
evolution  strategy  will  jump  ahe^d  as  the  result  of  a  large 


mutation  length 


2.3  Evolution  Strategy  Examples 

Several  examples  of  evolution  strategy  application  are  shown: ,  ;,*'2 


in  figures  (11)  and  (12)10’.11  As  a  simple  verification  test,  figurf  m 


(11B)  depicts  a  segmented  plate  having  5  degrees  of  freedom  that 
is  mounted  in  a  wind  tunnel.  The  objective  is  to  reduce  the  overall 
drag  of  the  plate  as  measured  by  the  momentum  defect  betv/een  the  " 
leading  and  the  trailing  edges.  Prom  the  initial  configuration,  g'; 
several  new  configurations,  or  cases,  were  created  by  random  mod¬ 
ifications  of  the  angle  of  each  segment  with  respect  to  an  adjacent 
segment.  After  1.40  cases,  the  segmented  plate  became  nenrly  planar, 
the  expected  result.  The  important  result  is  that  140  cases  rep-  • 
resents  only  4 .06*ifc0’"5$  of  the  total  possible,  345  million. 

The  second  example  is  an  attempt  to  modify  a  system  whose  op¬ 
timum  configuration  is  unknown^.  In  this  case,  shown  in  figure 
(12A),  the  objective  was  to  maximise  the  efficiency  of  a  two-phase 
supersonic  nozzle  constructed  from  a  number  of  discs.  The  optimal 
configuration  is  radically  different  from  conventional  designs  and 
verifies  the  ability  of  evolution  strategy  to  converge  to  unknown 
solutions.  It  is  also  interesting  to  note  that  flow  mixing  regions 
seen  in  the  figure  have  been  previously  suggested  for  improving  noz¬ 
zle  efficiency. 

The  third  example,  that  of  reducing  the  head  loss  for  a  90° 
bend  of  flexible  tubing,  also  has  an  unpredictable  result.  As 
depicted  in  figure  (12B)  the  final  share  is  a  subtle  change  from 
the  initial,  a  standard  circular  arc,  but  the  result  of  300  muta¬ 
tions  is  a  halving  of  the  head  lossll. 


V'ViSl 


2.4  Selection  of  Optimization  Method 

The  advantares  of  using  evolution  strategy  over  an  exhaustive 
evaluation  of  every  possible  configuration  are  clear,  except  perhaps 
for  the  case  of  a  system  having  only  one  degree  of  freedom.  The 
advante~es  of  evolution  strategy  over  a  ’strictly  determined 
mathematical  strategy’,  such  as  gradient  methods,  have  been  sum¬ 
marized  by  Rechenberg-'-O ; 

1  When  a  large  number  of  parameters  are  involved,  the 
evolution  strategy  attains  the  desired  result  more  rapidly 
than  the  more  familiar  strictly  determined  search  strategies, 
assuming  the  size  of  the  search  steps  to  be  the  same  in  both 
cases.  So  far,  this  could  only  be  proved  for  the  case  of  an 
n-aimensional  hyperplane  rising  in  any  arbitrary  direction. 

A  more  general  proof  is  being  attempted. 

Whereas  the  mathematical  search  strategies  used  so  far 


2 
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require  very  small  stops  (in  the  sense  of  the  t  rune  at  ion;  6  f '  ^thev 
Taylor  series  after  the  first  term),  the  evolution  method  can;> 
and  should  operate  also  with  much  larger  steps  which  exceed  the^ 
linear  region  in  the  neighborhood.  Talcing  larger  steps  sig-  •  •  ; 
nifies  r' 

(a)  in  many  situations  a  more  rapid  advance  towards  the 

desired  aim,  and  ■'-> 

(b)  a  shorter  time  to  decide  whether  the  step  taken  has  been  ^ 
successful  or  not.  (The  change  in  the  value  of  a  function  will 
generally  be  greater  in  the  case  of  a  large  parameter  change 
than  in  the  case  of  a  small  change.) 

3  A  so-called  "steady  signal”  in  the  measurement  of  the 
value  of  a  function  is  a  mathematical  fiction.  Disturbances  are 
always  present,  which  give  rise  to  errors  in  measurement.  The 
effectiveness  of  the  strictly  determined  mathematical  search 
strategies  is  markedly  reduced  by  small  errors  of  measurement. 

It  in  of  the  nature  of  the  evolution  method  (as  a  stochastic 
search  method)  that  small  random  errors  of  measurement 

cannot  have  a  decisive  influence  on  the  develooment  of  the 
process . 

4  Thorn  are  oases  in  which  the  more  familiar  mathematical 
search  strategies  must  seach  deadlock.  Such  cases  are  fre- 
cmently  observed  in  physics;  examples  are  hysteresis  phenomena 
and  locally  limited  extremal  values.  The  evolution  method  can 
generally  cone  with  such  situations  without  difficulty. 

5  The  algorithm  of  the  evolution  strategy  is  extraordin¬ 
arily  simple.  Thin  imnlios  that  the  effort  required  for  an 
autonisation  of  the  search  process  is  relatively  low. 

Since  the  work  associated  with  this  resort  was  strictly  numer¬ 
ical,  as  onnosed  to  experimental ,  the  small  error  effect  noted  above 
related  to  the  use  of  mathematical  strategies  could  not  have  any 
significance;  thus  an  argument  coxxld  still  be  made  in  favor  of  gra-'  " 
dient  methods,  Rochonberp,  however,  has  found  that  a  cross-over 
point  exists  at  which  search  effort  for  improved  configurations 
using  evolution  strategy  is  less  than  the  effort  required  for 
gradient  strategies?.  This  noint  is  for  systems  having  5  or  more 
degrees  of  freedom. 

The  system  of  interest  in  t^e  w^rk  described  here  i1-  a  A- 
element  airfoil  described  by  10  parameters  as  shown  in  figure  (13). 
Three  pivot  points  describe  the  relative  position  of  adjacent  ele¬ 


ments.  The  individual  geometry  of  an  element  is  given  with  respect 
to  a  local  coordinate  system,  with  the  origin  on  its  leading  edge. 
Pivot  points  1  and  2  are  fixe'1  with  respect  to  element  2  and  are 
variable  with  respect  to  elements  1  and  3.  Likewise,  pivot  point 


3  is  fixed  with  respect  to  element  3  and  variable  with  respect  to 
element  4.  The  two  degrees  of  freedom  for  each  of  the  pivot  points 
plus  the  relativp  deflection  betv/een  adjacent  elements  define  9 
degrees  of  freedom.  The  10th  degree  of  freedom  is  the  angle  of 
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attack  between  the  main  airfoil  a nd  the  fre^atream  velocity  vector* 
As  previously  stated,  the  objective  of  this  optimization  is 
to  obtain  the  best  envelope  of  Cl/Ca-vs-C^  values.  The  fitness  fac¬ 
tor  for  this  system  is  the  magnitude  of  C^/C^  for  each  value  of  C^, 
which  is  eq\iivalent  to  minimizing  for  fixed  Ci»  This  implies  a 
series  of  individual  optimizations  over  all  desired  values  of  0^. 

The  flexible  nature  of  evolution  strategy  allows  a  two-step  proce¬ 
dure,  however.  First  C^/C^  is  maximized  while  maintaining 
above  the  initial  value;  the  limit  of  this  maximization  is  a  point 
on  the  envelone.  The  second  steo  is  then  to  move  alonm  this  envoi- 
one  with  moderate  mutation  lengths,  selecting  those  configurations 
having  the  best  combination  of  0\/(3a  ano  Up .  it  in  because  of  this 
flexibility  couoled  with  the  ahovementioned  aspects  of  evolution 
strategy  that  this  method  was  chosen  over  gradient  strategies  an  a 
basis  for  multi-element  airfoil  optimization. 


3  CALCULATION  03?  AERODYNAMIC  COEFFICIENTS  .-v-'V  ;': 

3.1  Description  of  Program  Theory  and  Use  V.  ;;vMy . 

Aerodynamic  Coefficients  are  calculated  hy  th a  program'  ”  _  s 
detailed  in  reference  13.  This  program  is  composed  of  the  JoTlpwing 
task  areas:  .  :  - 

A.  Potential  Flow  Solution  '-dd 

B.  Ordinary  Boundary  Layer  Solution 

C.  Confluent  Boundary  Layer  Solution 

D.  Slot-Flow  Analysis 

E.  Combined  Solution 

The  inviscid,  potential  flow  solution  makes  use  of  the  distributed 
vortex  concept  with  the  vortex  singularity  comprising  the  funda¬ 
mental  solution  i  •>  the  Laplace  equation.  Airfoils,  arbitrarily  .s 
arranged  and  composed  of  from  one  to  four  elements,  have  surfaces 
approximated  by  a  closed  polygon  whose  linear  segments  are  represented 
by  distributed  singularities.  Airfoil  contours  are  limited  to  smooth, 
regular  shapes  with  sham  or  pointed  trailing  edges.  The  only 
restriction  on  the  slot  formed  by  two  adjacent  elements  is  that  the 
magnitude  of  flap  overlap  must  bo  greater  than  about  1/  of  the 
airfoil  chord. 

The  ordinary  boundary  layer  solution  is  comprised  of  mathe¬ 
matical  models  for  laminar,  transition,  and  turbulent  boundary 
layers  in  subsonic  flow.  The  laminar  boundary  layer  model  is  based 
on  the  basic  approach  of  Cohen  and  Reshotko  modified  to  suit  the 
needs  of  the  program.  Laminar  stall  criterion  developed  by  the  author 
of  the  program  predicts  the  formation  of  short  or  long  separation 
bubbles  or  bubble  burst,  a  flow  condition  which  causes  the  termination 
of  further  program  execution.  The  transition  model,  evolving  from 
the  instability  criteria  of  Schlicting  and  Ulrich  ,  establishes 
limiting  conditions  for  defining  the  position  of  transition  on  the 
airfoil.  Two  separate  mathematical  models  for  ordinary  turbulent 
boundary  layer  developmrnt  are  used.  The  first,  an  approximate 
model  developed  by  Goradia,  is  used  in  the  initial  iterative  calcu¬ 
lations.  The  second  and  more  accurate  model,  based  on  the  methods 
of  Nash,  determines  the  turbulent  boundary  layer  in  the  final,  vis¬ 
cous  solution. 


A  significant  feature  the  program  is  the  inclusion  of  a 
confluent  boundary  layer  model  that  reflects  ihe  merging  of  an  upper 
surface  boundary  layer  with  slot  efflux.  This  model,  developed  from 
the  experimental  and  analytical  work  of  Goradia,  accounts  for  the 


highly  complex  viscous  phenomena  associated  with  slotted  airfoils. 
Associated  with  the  confluent  boundary  layer,,  a  slot-flow  model 
is  defined  for  flow  between  slot  regions. 

The  final  viscous  solution  uses  an  iterative  technique  to 
combine  the  inviscid  solution  with  the  boundary  layer  calculations,. 
The  geometry  of  the  ’equivalent  airfoil’,  reflecting  local  boundary 
layer  thicknesses,  is  successively  defined  over  a  variable  number  of 
iterations  until  pressure  distributions  have  converged  to  a  steady 
condition.  As  the  work  described  in  this  reoort  was  more  of  a 
qualitative  nature  than  quantitative,  only  one  iteration  between 
viscous  and  inviscid  calculations  is  made,  in  the  interest  of 
reduced  computation  time.  A  comparison  of  predicted  and  experimental 
pressure  coefficients  using  the  program  for  a  two-element  airfoil 
appears  in  figure  (15)*  A  listing  and  description  of  program  input 
is  given  in  appendix  (i)  for  the  4-element  airfoil  considered  in 
this  work. 
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4  OPTIMIZATION  PROCEDURE 
4 .1  Geometry  Modifications 

The  choice  of  evolution  strategy  reouired  the  preliminary 
development  of  the  following: 

A.  Scheme  for  Random  Geometry  Modifications 

B.  Mutation  Length  Control 

C.  Criterion  for  Convergence  to  Optimum 

D.  Initial  Airfoil  Geometry 

Of  fundamental  importance  to  the  effectiveness  of  evolution  strategy 
is  the  combination  of  large  and  small  mutation  lengths.  If  only 
small  siaes  are  used,  the  number  of  improved  cases  will  be  high,  but 
the  rate  of  progress  will  be  slow.  On  the  other  hand,  if  only  large 
sizes  are  used,  the  optimum  may  be  entirely  skinned  over.  For  this 
reason,  new  cases  were  divided  into  3  groups,  one  having  all  mutation 
lengths  of  a  base  value,  the  second  group  having  lengths  of  50 fc  of  the 
base  value,  and  the  third  groun  having  lengths  15055  of  the  base., 

The  control  of  the  value  of  this  base  mutation  length  will  be 
described  below,  For  each  new  geometry,  then,  18  new  configurations 
are  generated,  divided  into  3  groups  of  6  cases  each.  The  choice  of 
18  cases  per  set  was  somewhat  arbitrary,  but  it  was  felt  that  this 
would  be  a  good  trade-off  between  effort  and  resulting  improvement. 

As  the  optimization  procedure  was  carried  out,  it  was  found  that  18 
cases  was  a  sa-fe  choice,  as  the  program  used  to  determine  Ci  and  C^ 
had  certain  limitations.  In  experiment,  a  smaller  number  of  cases 
per  set  might  be  desiroablc. 

All  mutations  wore  generated  by  a  computer  program,  whose 
listing  appears  in  appendix  (ii)  and  is  described  by  its  logic 
diagram  in  figure  (ISA).  The  random  choice  required  by  this  program 
is  made  using  a  random  nufnber  generator  that  determines,  for  each 
case  which  of  the  10  parameters  is  to  be  modified,  and  in  the  event  of 
a  modification,  whether  the  increment  is  positive  or  negative.  As 
seen  in  figure  (16C),  two  columns  of  mutation  lengths  are  printed. 

The  first  gives  the  three  lengths  used  to  modify  the  x,y  position 
of  pivot  points.  These  values  represent  a  percent eye  of  the  range 
defined  for  each  pivot  point,  as  shown  in  fipire  (163).  The  second 
column  is  degrees  of  mutation  for  flap  deflection  and  angle  of  attack. 

The  Random  number  generator  is  a  computer  system  function  that 
returns  an  even  distribution  of  real  values  in  the  range  0. 0-1.0. 

A  description  of  its  algorithm  is  given  in  appendix  (iii). 
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4.2  Mutation  Length  Control  * 

As  described  above,  three  modification  lengths  are  used'  for  ' 
each  new  set  of  configurations.  The  control  of  the  base,  or 
middle  length  is  critical,  as  rapid  convergence  to  an  optimum  requires 
an  expansion  of  mutation  size  far  from  the  optimum  and  a  compression 
close  to  it.  In  this  optimization  problem,  the  base  length  for  a 
new  set  was  controlled  by  the  best  case  from  the  previous  set. 

For  instance,  if  case  1  ,  a  member  of  the  third  mutation  length  group, 
is  the  best  in  its  set,  its  mutation  length  will  be  used  as  the  base 
value  for  the  next  set.  In  this  way,  increasing  lengths  will  be 
favored.  On  the  other  hand,  if  case  3  is  the  most  improved,  smaller 
lengths  will  be  favored  in  the  next  set.  Using  this  scheme,  the 
optimisation  strategy  self-adjusts  to  the  distance  from  convergence, 

4.3  Convergence  Criterion 

Concerning  the  auestion  as  to  when  convergence  has  been 
reached,  a  criterion  is  used  that  recognizes  that  mutation  lengths 
will  automatically  compress  as  the  point  of  convergence  is  approached. 
When  the  modification  of  an  improved  configuration  fails  to  increase 
Ci/Cfl  for  a  minimum  Cx,  the  airfoil  has  ’ approximately  converged’,. 

For  absolute  convergence,  an  additional  sot  of  configurations  is  gen¬ 
erated,  based  on  the  best  case  of  the  set  that  failed  to  improve 
over  the  point  of  ap-nroximate  convergence.  If  avain  there  is  no 
improvement,  convergence  is  said  to  exist. 

The  procedure  proposed  to  determine  the  envelope  of  aero¬ 
dynamic  coefficients  is  comprised  of  the  following  two  basic  steps: 

1,  Maximize  Gn/C^  while  maintaining  Cx  greater  than  Cxmin» 

the  value  obtained  from  the  initial  geometry 

2 .  Move  across  the  envelope  from  the  maxima  point  found  from 

step  1  to  the  maximum  Cx  limit,  when  flow  separation  occurs 

Results  from  this  procedure  are  presented  in  section  5, 

4.4  Initial  Airfoil  Geometry 

The  initial  airfoil  configuration  was  an  arbitrary  but 
realistic  placement  of  elements  so  as  to  guarantee  good  flow  condi¬ 
tions  through  the  slots.  The  value  of  the  10th  parameter,  angle  of 
attack,  was  determined  from  the  polar  shown  in  figure  (14).  .ftn 
angle  of  0.5°  was  used  as  it  resulted  in  the  maximum  Cx/C^.  This 
initial  configuration  and  angle  of  attack  thus  set  the  value  of 
Ci  •  at  2.“1.  In  all  subsequent  program  calculations,  mach  number 
and  Reynolds  number,  referenced  to.  an 'airfoil  chord  of  350mm, 
were  0.125  md  1.26*10^,  respectiv°ly.  Element  profiles  were 
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taken  from  reference  14,  with  runr 
the  reouirement  of  smooth  contours 


5.1  First  Three  Configuration  Sets 


The  first  step  of  the  optimization  procedure,  the  maxim-  '  1 

ization  of  Ct/C.  with  the  constraint  of  Ci  .  ,  was  started  at 

J-  a  •'•mm  % 

the  initial  geometry.  Figure  (17)  shows  the  scatter  of  data.)  * 
points  for  the  first  three  sets  (54  cases)  of  airfoil  configurations,  \ 
Some  points  have  been  omitted  for  clarity.  Circled  data  points  .are  j 
’best  cases  for  each  set;  these  cases  serve  as  the  bases  for  subse-  "*’<  j 
‘quent  modifications.  ’  The  automatic  mutation  length  control  reacts 
well  in  these  first  sets,  far  from  the  point  of  convergence,  with  '*v  j 

an  increase  in  base  lengths  between  sets  one  and  two  and  a  constant  \ 

base  length  between  sets  two  and  three.  As  exnected,  points  scattered  ; 
in  the  immediate  vicinity  of  the  initializing  case  wore  from  the  first 

j 

mutation  length  group,  and  points  scattered  further  away  were  from  the  ! 

second  and  third  groups.  As  seen  in  figure  (17),  the  ratio  Cx/C^  .  | 

* 

apnoars  more  sensitive  to  airfoil  modification  than  Cx.  This  effect  j 
continued  for  all  other  sets,  with  most  points  located  in  a  narrow  : 

Cx  band  between  2.0  and  3.0  and  a  scatter  of  other  points  in  a  rel¬ 
atively  low  Cx/Cd  -  low  Cx  region.  The  tabulation  of  all  configurations 

(37  sets,  656  cases)  is  given  in  Appendix  (iv)  together  with  initial  ! 
and  final  configuration  data.  V’*f  i  j 


,5.2  Description  of  Complete  Optimization  .  -  .  --  • 

Figure  (18)  shows  the  result  of  continued  maximization  of 
Cl/Cfi,  plotting  tho  best  case  for  each  of  37  sets.  Up  to’  set  11,  ... 
mutation  lengths  generally  expanded  or  remained  constant.  Sets  - 
12-14,  however,  resulted  in  a  compression  of  lengths,  ana  convergence 
was  thought  to  be  imminent.  Beyond  set  14,  lengths  began  to  expand 
again  until  the  base  value  stabilized  at  approximately  2e/'  and  1° 
for  pivot  point  location  and  deflection  angles,  respectively.  This 
effect  was  considered  to  be  a  'small  barrier'  in  the  optimization  . 
process,  which  the  evolution  strategy  was  able  to  step  over  with  the 
larger  of  the  three  mutation  lengths.  A  gradient  strategy-based 
technioue  would  have  converged  to  such  a  barrier,  or  local  maximum. 

Beyond  set  22,  the  accuracy  of  tho  aerodynamic  coefficients 
became  questionable,  resulting  in  erratic  behavior  of  the  optimiza¬ 
tion  process.  When  the  process  was  continued  without  the  Cxm^n 
constraint,  large  changes  in  Cx/C(j  and  Cx  occurred  for  relatively 
small  geometry  modifications.  As  an  attempt  to  return  to  smaller 
variations  of  coefficients,  the  process  was  restarted  at  set  24 


i 

i 

< 


with  a  forced  reduction  of  step  sizes,  as-  shown  by  the  dotted  lines1 
of  figure  (18).  This  reulted  in  sets  27  and  28  which  again  were 
characterized  by  large  variations  of  coefficients  for  small”  modifi¬ 
cations.  This  sequence  of  forced  length  compression  and  optimiza¬ 
tion  restart  was  continued  to  the  37th,  but  by  this  point  values 
of  C-j/C^  were  unrealistically  high,  even  for  the  simplified  two- 
dimensional  airfoil  model  without  induced  drag  or  body  interference 
effects.  Indeed,  during  the  process  of  optimization,  numerous 
configurations  resulted  in  unsuccessful  program  termination  or 
extremely  high  values  of  Cq/Gg.and  in  several  cases,  negative  C&! 

The  reasons  for  this  inaccuracy  are  thought  to  be  one  or  a  com¬ 
bination  of  the  following: 

1  Invalid  geometry;  i.e.,  overlap  less  than  lje  of  chord 

2  Confluent  boundary  layer  or  slot  flow  model  failures  for 
small  slot  cross-sections 

3  Inadeounte  verification  of  program  by  authors  for  4-element 
airfoil  configurations  (note-the  freciuoncy  and  severity  of 
inaccuracies  were  substantially  les<-  for  the  preliminary 
optimization  of  a  2-element  airfoil) 

In  the  further  presentation  of  results,  rets  21-22  are  taken  to  be  the 
apnroximate  limit  of  data  accuracy,  although  all  points  ore  shown. 

Gl/Cd  Gl  arG  shown  with  respect  to  set  number  in  figure 
(19).  The  points  again  are  from  the  best  cases  for  each  set.  As 
seen  in  this  figure,  the  result  of  22  sets  of  configurations  and 
296  total  cases,  C^/Cfl  has  been  increased  from  4.28  to  36.6,  a  factor 
of  8.515.  At  the  same  time,  Cq  has  been  increased  from  2.31  to  2 . 5 9 » 
a  result  that  evolved  by  favoring  the  configuration  having  the  larger 
Cl  when  two  or  more  had  aonrocimately  the  same  Cj/Cq.  Increasing 
Ct,  however,  was  not  the  primary  objective  of  the  optimization. 

Because  of  the  limit  of  data  accuracy,  the  optimization  process 
was  terminated  at  the  37th  set,  thus  making  it  impossible  to  find 
the  envelope  of  aerodynamic  coefficients  and  thus  carry  out  the 
second  ste^  of  the  proposed  optimization. 


5.3  Initial-Final  Configuration  Comparison 

?->  "ure  (20)  shows  polars  generated  for  the  initial  configur¬ 
ation  and  the  configurations  resulting  from  the  process  of  Ci/C<j 
maximization.  The  resulting  dramatic  improvement  in  aerodynamic 
performance  is  apparent.  During  the  process  of  ootimization,  a  num¬ 
ber  of  configurations  produced  values  of  Ci/Ca  and  Ci  that  when  plot¬ 
ted  fall  to  the  right  of  the  polar  for  maximized  Ci/C^.  A  polar 
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was  generated  from  tlie  point  that  was  furthest,  to  the  right  and  is; 
denoted  as  the  ’optimized  configuration’.  This  configuration,  it 
is  interesting  to  note,  was  a  member  of  the  optimum  •%-/<%•  'set ,; 
number  22.  The  dotted  lines  have  been  added  to  these  polar s  to  show 
the  real  flow  behavior  for  multi-element  airfoils  such  as  shown  in 
figure  (2). 

A  comparison  of  initial  and  final  configurations  is  given  in 
figure  (21).  A  noticeable  effect  of  the  optimization  is  the 
reduction  of  slot  area  and  improved  contour  smoothness  of  the  trailr 
ing  edge  flaps.  Both  of  these  effects  have  been  shown  experiment ally 
to  improve  aerodynamic  performance.  The  same  changes  did  not  occur 
for  tho  leading  edge  flap.  It  was  found  that  the  multi-element  .  ► 
airfoil  program  was  very  sensitive  to  reduction  of  slot  area  of  the 
leading  edge  device  when  both  trailing  edge  flaps  were  present.  It 
is  suspected  that  the  reasons  for  this  sensitivity  arc  the  same  as 
those  mentioned  above. 

A  comparison  between  tho  optimum  C]/Ca  configuration  found 
using  evolution  strategy  and  that  suggested  by  the  empirical  reco¬ 
mmendations  of  reference  5  was  attempted,  but  tho  slot  geometries 
required  were  not  capable  of  being  modelled  by  the  airfoil  program. 
Any  change  from  these  recommendations  so  as  to  suit  program  input 
requirements  would  have  resulted  in  a  meaningless  comparison. 


CONCLUSIONS- 

*  ( 

*  ' :  ■  , 

Although  the  inaccuracy  of  the  program  used  to  determine  .aero- 
•dynamic  coefficients  precluded  the  full  solution  of  the  airfoil  per¬ 
formance  envelope,  the  problem  of  maximizing  C^/C^  was  easily  handled1 
by  the  evolution  strategy-based  optimization  technique.*  A  review  of 
existing  techniques  showed  that  evolution  strategy  is  the  best  choice 
for  systems  described  by  5  or  more  degrees  of  freedom,  due  to  its 
relative  incomplexity,  flexibility  to  handle  a  wide  variety  of  prob¬ 
lems,  and  to  the  knowledge  that  the  converged  solution  is  an  absolute 
optimum  within  the  range  of  variation  of  these  degrees  of  freedom. 

The  problem  experienced  with  the  use  of  the  above-mentioned 
program  reinforces  the  opinion  that  as  of  yet  numerical  solution 
methods  are  severely  limited  in  their  range  of  application. to  complex 
engineering  systems,  such  as  for  a  multi-element  airfoil  with  its 
large  number  of  configurations.  The  value  of  numerical  methods, 
however,  lies  in  their  ability  to  provide  preliminary  results  to  be 
used  as  a  starting  point  for  further  design  refinement  in  the 
laboratory.  Of  course,  as  long  as  the  limits  of  numerical  methods 
are  not  exceeded,  experimental  optimization  may  bo  well  predicted. 

The  next  logical  step  related  to  the  work  described  here  is  the 
application  of  evolution  strategy  to  optimisation  in  a  wind  tunnel. 

A  setup  is  envisioned,  similar  to  those  of  references  7  end  8,  where 
modifications  to  the  airfoil,  analysis  of  aerodynamic  coefficients, 
and  optimization  procedure  are  carried  out  by  on-line  computers. 

The  importance  of  using  the  wind  tunnel  is  that  all  configurations  of 
an  airfoil  are  capable  of  being  tested,  as  opposed  to  computer  model- 
imposed  restrictions  such  as  slot  geometry  and  separation-free  flow. 
An  important  note  is  that  multi-element  airfoils  often  are  designed 
to  perform  with  separated  flow  and  negative  overlap,  two  conditions 
that  can  not  exist  for  the  successful  use  of  the  program  of  refer¬ 
ence  13. 

The  effect  of  increasing  Ci/C^  through  the  modification  of 
slot  geometries  while  at  the  same  time  maintaining  adeouate  Ci  is 
best  seen  for  the  design  trade-off  between  take-off  and  cruise  flight 
Increased  C3/C3  at  take-off  allows  increased  rate  or  angle  of  climb, 
decreased  overall  distance  to  reach  a  required  screen  height, 
increased  take-off  weight  for  the  same  runway  requirement,  or  a 
combination  of  the  three.  An  increase  in  Ci/C $  at  take-off  also  has 
a  stronm  influence  on  cruise  flight  efficiency  througn  the  reduction 
of  the  installed  power  needed  for  take-off. 
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By  using  an  evolution-based  strategy.,  -the  converged  optimum 
is  an  ’absolute*  optimum  within  the  range  of  variations  of  the 
degrees  of  freedom  of  a  particular  system.  The  question  of  whether 
convergence  has  been  obtained  or  not  is  thus  eliminated;  this  was 
not  the  case  for  the  design  of  the  Boeing  737.  This  is  seen  i'n 
figure  (4)  by  the  3  increments  in  Ci  for  its  high-lift  devices i 

iuwa 

At  each  point,  the  design  was  probably  considered  to  be  ah  optimum! 

Other  areas  of  aerodynamic  design  that  should  lend  themselves 
well  to  evolution-strategy  optimization  are  the  wing-fuselage 
.junction,  empennage,  supercritical  airfoils,  engine  nacelles,  and 
fuselage  rear  up-sweep.  These  design  areas  are  currently  approached 
using  procedures  similar  to  those  of  flap  design,  the  wind  tunnel 
testing  of  a  relatively  small  number  of  different  configurations. 

Each  of  these  problems  is  described  by  many  degrees  of  freedom,  a 
characteristic  that  has  been  shown  to  be  well-suited  to  evolution 
strategy.  The  technioue  described  in  this  work  is  by  no  means  restr¬ 
icted  to  aerodynamic  design,  however.  A  wido  range  of  applications 
can  be  imagined;  the  examples  included  in  this  report  are  only  a  few 
preliminary  cases  whoso  results  show  promise  for  this  new  technique  to 
be  used  as  a  powerful  tool  in  the  design  process. 

One  can  consider  the  development  of  evolution  strategy  as  a 
result  of  man  learning  from  his  observations  of  nature.  Figure  (22) 
depicts  what  might  result  if  nature  learns  from  the  aerodynamic 
achievements  of  man. 


19 


REFERENCES 

1.  Wimpress,J,:  "Shortening  the  bakeoff  and  Landing  Distances  of 
High  Speed  Aircraft,"  the  Boeing  Company,  Ren-o on,,  Washington, 
D6-16168,  June  1965. 

2.  Steiner,  J.,  et.al.,  "Case  Study  in  Aircraft  Design:  The  Boeing 
727/'  AIAA  Professional  Study  Series,  Sept.  1978. 

3.  Bruner,  G.,  "Le  Boeing  747,"  DOC-AIR-ESPACE ,  No. 112,  Sept.  1968, 
PS.  7. 

4.  Technical  Staff,  AEDC  (USAF)  and  ARO,  Inc.,  edited  by 
J.D.  Whitfield,  J.P.  Hartin,  and  S.R.  Pate,  "Aerodynamio 
Testing-a  Look  at  Future  Requirements,"  AIAA  No. 78-765,  1978. 

5.  L.iungstrbm,  B.L.G.,  "Experimental  High  Lift  Optimisation  of 
Multiple  Element  Airfoils,"  Proceedings  of  the  AGARD  Conference 
on  V/STOL  Aerodynamics ,  AGARD-CP-143 ,  Apr.  1974,  ?o.  13-1  -  13-16 

6.  Rettie,  I.K.,  "Air  Transportation  for  the  1980’s:  Aerodynamic 
Design,"  from  collected  papers  of  the  1975-76  Seminar  Series 
supported  by  the  Minta  Martin  Fund  for  Aeronautical  Research, 
University  of  Maryland,  College  Park,  Department  of  Aerospace 
Engineering,  Juno  1976,  pp.  115-152. 

7.  Levinsky,  E, ,  Schappellc,  R.,  and  Pountney,  S.,  "Airfoil 
Optimization  Through  the  Adaptive  Control  of  Camber  and 
Thickness ,  Phase  II:  Transonic  Wind  Tunnel  Tost  and  Program 
Summary,"  General  Dynamics  Corporation,  Convair  Division, 
CASD-NSC-75-004 ,  Sept.  1975. 

8.  Levinsky,  E.S.,  Palko,  R.L.,  "Semispan  Wind  Tunnel  Test  of  a 
Computer  Controlled  Self-Optimizing  Flexible  Technology  Wing," 
AIAA  paper  78-786,  1978. 

9.  Hicks,  K. ,  Kurman,  E. ,  and  Vanderplaats ,  G.,  "An  Assessment  of 
Airfoil  Design  by  Numerical  Optimization,"  NASA  TM  X-3092 , 

July  1974. 

10.  Reehenberg,  I.,  "Cybernetic  Solution  Path  of  an  Experimental 
Problem,"  Royal  Aircraft  Establishment  Library  Translation 
No.  1122,  Aua.  1965. 


-  20  - 

REFERENCES  (cont’d)  .  ;V 

11.  Rechenberg,!.,  "Evolutionsstrategie,  Optimierung  Technischer 
Systeme  nach  Prinzipien  der  Biologischen  Evolution,”  No.,  15  t 
of  Problemata  Series,  Friedrich  Fromraann  Verlag,  Stuttgart- 
Bad  Cannstatt,  1973. 

12.  Discussion  between  Professors  i.  Rechenberg  and  R.  Stuff, 
Technischen  Universitat  Berlin,  Feb.  1980. 

13.  Stevens,  W. ,  Goradia,  S.,  and  Braden,  J.,  "Mathematical  Model 
for  Two-Dimensional  Multi-Component  Airfoils  in  Viscous  Flow,” 
T.ockheed-Georgia  Company,  Marietta,  Georgia,  NASA  CR-1843  and 
Supplement  to  NASA  CR-1843 ,  July  1971. 

14.  Hertha,  K.  and  Scheerer,  J.,  ”Konotruktion  des  ZKP-Panelhalbmodell: 
Kesserschmitt-3blkow-31ohm  Gmbh,  Hamburger  Flugzeugbau,  LFK  7512, 
Ergebnisbericht  Nr.  19,  Arrbeitspaket  :E-I.4 ,  April  1977.. 

15.  Schwefel,  Hans -Paul ,  "Numerische  Optimierung  von  Computer- 
Kodellon  mittels  der  Evolutionsstrategie,  volume  26, 
Interdisciplinary  Systems  Research  (ISR  26),  Birkh&uscr  Verlag, 
Basel,  1977. 


*  *  '»i*Tr»*r^3^*.''  WWn*  -9 


FLAPS  UP 


O  727  TAKEOFF 
&  727  LANDING 


THEORETICAL  DRAG 
FLAPS  UP 

— t-, APPROXIMATE 
\  ‘  •  ENVELOPE 


40  FLAPS 


0.0  0.5  10  15  2.0  2.5  3.0 

LIFT  COEFFICIENT- CL 


BOEING-727  DRAG  POLARS 


FIGURE  2 


MULTI-ELEMENT  AIRFOILS 


FIGURE  3 


DEFLECTION 


L.E.  SLAT  GAP  &  DEFLECTION  EFFECTS 


T.E.FLAP  GAF  &  OVERLAP  EFFECTS 


COMBINED  SLAT  AND 
FLAP  EFFECTS 


T 


FOUR  ELEMENT  AIRFOIL  OPTIMIZA 
TION  BY  EMPIRICAL  METHOD 


FIGURE  5 


10  w 


F* - 100  YEARS 


SST 


SHUTTLE 


— 1 .01 
2000 


TIME,  YEAR 


WHAT  PRICE  ACCURACY  ? 

•  •  •  •  W10  CONFIGURATIONS 

■  ■  ■  ■  1  MEASUREMENT/ 30  SECONDS 

■  ■  ■  TOTAL  TIME  9.51 -10s  YEARS 

■  ■  ■  TOTAL  COST  $B.33-1010 


WIND  TUNNEL  TEST  PROGRAM 
TRENDS 


FIGURE  6 


THREE-DIMENSIONAL  FLEXIBLE 
AIRFOIL  OPTIMIZATION 


FIGURE  8 


-1.0 


-.a 


0  10  20  30  40  50  60  70  80  90  100 

X/C,  PERCENT  CHORD 
—  INITIAL  AIRFOIL ,  Crf  0.00*3,  V*0.721 
— FINAL  AIRFOIL ,  CqbO.0009,  \ 1=0.756 
VOLUME  MAXIMIZATION ,  0^*0.001,0^0.0, 


X/C,  PERCENT  CHORD 


—INITIAL  AIRFOIL ,  CD=0.0524 ,  \ 1=0.660 
—FINAL  AIRFOIL ,  CD=0.0007,  V=0.654 
DRAG  MINIMIZATION ,  V  20. 65,0^0.0,  M=0.d5 


TWO-DIMENSIONAL  TRANSONIC 
AIRFOIL  OPTIMIZATION,  GRADIENT  ST. 


FIGURE  9 


t/c,%  CHORD  x  i/c,  %  CHORD 


AIRFOIL  WITH  10  PARAMETERS 


FIGURE  13 


POLAR  FOR  INITIAL  GEOMETRY 


FIGURE  74 


COMPARISON  OF  EXPERIMENTAL  AND  PREDICTED  PRESSURE 
DISTRIBUTIONS  FOR  NACA  23012  AIRFOIL  WITH  T.EDGE  FLAP 


parameter 
=  initial 


call 

randu 


call 

randu 


add  Jo 
init .  aeom . 


4.  PROGRAM  FLOWCHART 


LENGTH  GROUPS — 

10  PARAMETERS  - 
INITIAL  GEOMETRY ■ 


LENGTH  GROUP  NUMs 

CASE  NUMBERS 
PAR3 -~y  /• 


fM  .<ui'lr  ICA  f  lo  4  O - -S£T  NUMBER 

1 1  *  I'll'-il')  IHi  Ntf.Ci  1  sflSts  6 


%i£y- 


B .  WE  IV  OF 
PARAMETERS 


1M.H  NO- 

i  i  1 1  ** 

uiiii'i  r 

1 

1  ,  •> N 

4*,6  • 

4 

y.  »  ^ 

1  .  <  > 

i 

v.  16 

2 , 0  v 

MM  PAM  i  AK6  I  Adb  P Ah  7hAK  o  PAM  P*hlU 
l.»'A^  J  l  l  I  1  1  1  1  |  1 

-O.VO-K  ,v*>  0,<*»  O ,  C>u  16.mO  O.to  O,*©  l‘j.00  O.’jO 

J  ’•  v-r  ,*  C.*-M  A  ?  > 

t  l  1.60  -0 ,  -.S-  2  * ,  bo  o./6  4*,<>0  •  *>,yc4  '**,,‘3  O.vO 

cm  j  .'v  Ci*s  <  u/u;  *',*  66  Cm  ~ 

l  2  1.1*0  46*  JO,  )•>  o.W  o.ts  M .  6o  o,0v  0.J3  16,410  0,t>o 

Cu  VwJ  CM  C>  v!.'  CL/vP*  /  c**  — P  6/4 

I  1  I.n©  -0.66-10,.  )  v , *•  •>  M,6f>  G.w  o,]3  16.  "0  u,6o 

vl-  .  ;  h  '*'■-£  5: 0  cis  -  ,  7*0 

flj  1  l.M  •>(  •< ,  4  |  »',«.*  I’J.W.1  •  o  ,  t  ©  li»*0>7  0.60 

v-r  Cl  -  '  VI  •  ..\‘U2  UKUS  /,  7  J.  l>*  -f. *1*1*4 


,  S<>«  . 

2  *V  CM 

•U,4l*  U>  , 
»  - 5-  Cl  = 

1 1  , 


-u.So-  lo, 
?, ci-s 

-0,  «»0-  1  I  , 

,0  C“- , 

-M.t  <#•  11 . 
■  v  n  C  = 

•'•.S'1-?  4  . 

r  ;  c  1  - 


■  j 4  k^:iac 

.  <  O 

/w  Cl'*  *  5  ^  > 


l  .  /  >  9  j  m  >S'  O. 

ui  '  £7 

i . « I  -J'4  >6-. .  _>i, 

»**•**-'  v  -r- £$  . 

I  .  2  v< .  <1,  u  , 

C  L  '  ►  -  •*  V*:  /  ‘ 


.©*>  14,60  0 

=  <;  7?  9  C  ’♦  * 

.**6  16,60  -o 

=  4.  /  V  C«s  ' 

,  M>  1  6 , 1  'j  0 

s  r  76  C»* 

.00  1  J  ,  v  O  0 

!  7  *i  C*  ; 


^.250  U./U*  7  C*  = 

j  •*  ,i*o  0 

•  •  6/  v.L/Ck  s  6  w  ^  V  *•  = 

,  i<  *  Mu  3,oo  l  •> ,  uii  0 
S  ASS  VL/V  I  *;!•!•{  c 

p1'*1  »  .  ?%  »,VU  ( *>.  1»0  1* 

^  s  ?i  ii/i'vs  j  a  os 


,u  I V .  k  0  l* 

=  4  0/  0  = 


It.vu  <* 

;  ?.>*/  C*s 

,*>J  l  l.uii  »u 

-  7J  V-s 

,  t»u  1  j  ,  / )  v> 

-  J  *  J.  C  -  s 

,0"  I  J  , V  w 

-  J  r  >  v-  = 

.  ©  »  1  J  ,  1  o  <* 

-2  *  a  c ^  s 

.00  I  *J.»*W  -  . 

=  <r  c  *  = 


,0d  O , 2 0 

'•/,  /3 

,09 

-/.^O 

,  1  b  0  t  t  n 

-/>  c« 

,!»*)  0,2  0 

-  0.Z&/ 

.1*  *>,:© 

-/ / f 

,  **o  0,1© 
-O  b  7.*i 

.00  »‘,^o 

-*/  22 

«-,i© 

-  /  /  i 

, '  2 

-;;-sy 

,  1‘0  0,^0 

- 

,  ?  2  0  ,  ^  Cl 

-/  7  2 

.<.*  **.4» 

. ’J./rt 

-A  24 


»  ?*5  1*>.C0  o.v 

s  s:  i-s  - 


16,00  1,00 
16,00  O.OV 
14,00  *1,60 

46,00  0,60 

14,00  J,60 

16.00  0,6y 

l**.  00  o ,  60 
16,00  1.60 
16,00  0,6-* 
U.t;o  0,60 
13.00  g,‘0 

IS.vO  **.6y 
4  S ,  y  0  O.60 

li.lO  2.6o 


C.  PROGRAM  OUTPUT 


GEOMETRY  MODIFICATION  PROGRAM 


FIGURE  16 


o  o 


OPTIMIZATION  PROCESS  FOR 
FIRST  THREE  SETS 


FIGURE  17 


100.0 


FIGURE  19 


C,  ANdB-VS.  set  number 
L  CD 


COMPARISON  OF  INITIAL  AND  FINAL  CONFIGURATIONS  I  FIGURE  27 


EVOLUTION  STRATEGY  -  THE  SKY  IS  THE  LIMIT ....  I  FIGURE  22 


A1  - 


^EI'TDIX  (i)  A?RODv?T/-j',.IC  CQRF^CTPiX'1  PRC/g/g  THPUg 


100 

••BB-PROPILE 

Al-725-1 46 

SLOTTED 

KRUEGER 

DOUBLE-SLOTTED  TE  FLAP  OPTIMIZATION 

Oi  0  0 

4  85 

1 

J 

300 

400 

1  23 

1.6100 

23 

-0.4200 

500 

0.0000 

0.0344 

0,0689 

0.1378 

0.2067 

0.2756 

0.3445 

0.4134  • 

bOO 

0.5512 

0,bb90 

0.8268 

0.96-16 

1.1024 

1.2402 

1.3780 

1.51S7 

7u0 

1.0535 

1.7913 

1.9291 

2.0670 

2.2047 

2.3425 

2.4459 

800 

0.4134 

0.5994 

0.6683 

0.7441 

0,7992 

0.8543 

0.8819 

0.8957  i 

900 

0.9232 

0.9232 

0.9163 

0.9026 

0.B888 

0 , 8  u  1 2 

0.8406 

0.8130  ; 

1000 

0,7923 

0.7579 

0.7303 

0,6890 

0.b407 

0,5925 

0.5512 

. 

1100 

0.0000 

0,0344 

0,0689 

0,1378 

0,2067 

0,275b 

0.3445 

0,4134  . 

1200 

0.5517 

0.6890 

0 , 8  2  b  8 

0.964b 

1.1024 

1.2402 

1.3780 

1.5157 

1300 

1  ,b535 

1.7913 

1.9291 

2.0670 

2.2047 

2.3425 

2.4459 

1 

1400 

0.4134 

0,2070 

0. 137H 

0,0689 

0.0207 

0.0069 

0.0000 

0.0207 

1500 

0 , 09o5 

0.2070 

0.2894 

0.3445 

0.3U58 

0.4272 

0, 4b 1 6 

0.4961  . 

1 600 

0.5236 

0.5443 

0.5512 

0.5512 

0.5512 

0.5512 

0.5512 

t 

1700 

2  37 

3  7 

1600 

-0.5000 

0.0000 

12,0000 

0.0000 

1000 

0.0000 

0.0J44 

0.1033 

0.1722 

0.2411 

0.3100 

0.3789 

0,4478 

2000 

0 , 5 1 b7 

0.565b 

0.6614 

0.8364 

1,0486 

1.3049 

1.6150 

1.9898 

2100 

2.4145 

2.9929 

3.6571 

4,3543 

5.3986 

6,4433 

7.1406 

7.7358 

2200 

8.2161 

8. 6122 

8,9264 

9,1951 

9 ,  b  8  7  0 

10, 1004 

10.6033 

11.2110 

2300 

1 1,5596 

11.9110 

12.2252 

12,5091 

12.5394 

2400 

0.0000 

0.1654 

0 , 2b  3  2 

0,3236 

0.3707 

0.4079 

0.4396 

0.4685 

2500 

0.4919 

0.5154 

0.5374 

0.5843 

0.6339 

0 . 6b  48 

0.7372 

0,789b 

2bOO 

0.8433 

0.8929 

0,9384 

0,9701 

0.9894 

0.9  756 

0.9467 

0,90e7  , 

2700 

0.H684 

0.8254 

0.7896 

0.7551 

0. b862 

0.6228 

0.5374 

0.4258 

20  00 

0.3596 

0.2921 

0,2301 

0.1791 

0.1722 

0.4478 

2900 

0.0000 

0,0344 

0,1033 

0.1722 

0.2411 

0.3100 

0.3789 

3000 

0.5167 

0 , 5b5o 

0 ,60  1  4 

0,8364 

1,0486 

1.3049 

l.blSO 

1.9898 

3100 

2.4445 

2.9929 

3.0571 

4.3543 

5.3988 

6.4433 

7,1406 

7.7356 

3200 

8.2181 

8.6122 

8.9264 

9,1951 

9 ,  b8  70 

10.1004 

10,b033 

11.2110 

3300 

1 1 ,559b 

11.9110 

12.2252 

12.5091 

12.5394 

-0.5434 

-0.5787  i 

3400 

0.0000 

-0.2267 

-0.3514 

-0.4203 

-0.4713 

-0.5140 

3SuO 

-0.6063 

-0.631  1 

-0.6559 

-0,7069 

-0.760b 

-0,8157 

-0.8681 

-0,9163 

JoOO 

-0.9577 

-0.9894 

-1  .0031 

-0,9921 

-0.9398 

-0.6392 

-0,246.' 

-0.7427 

-0.6476  ! 

37  00 

-0.5650 

-0.4974 

-0.4437 

-0.3982 

-0.3142 

-0.1695 

-0.0897  t 

3600 

-0,0551 

0.0344 

0.1378 

0, lb54 

0.1722 

1 

3900 

2  2  3 

23 

-0,1376 

1 

4000 

-0.1000 

0.2200 

2  •  06b9 

0.2136 

0.4134  ’ 

4100 

0.0000 

0.0138 

0.0413 

0.0896 

0.1516 

0.275b 

4200 

0.5512 

0,6890 

0.8268 

0.9646 

1.1024 

1.2402 

1.3780 

2,4803 

1.5157  t 

4300 

1.6535 

1.7913 

1.9291 

2 , 0669 

2.2047 

2.3425 

] 

4400 

0.0000 

0 , 09b5 

0.1376 

0,1929 

0.2343 

0.275b 

0.3031 

0.3445 

4500 

0.3720 

0.3858 

0.3927 

0,3856 

0,3/89 

0.3563 

0.3307 

0.3031  j 

4b00 

0,2756 

0,2343 

0.1929 

0,151b 

0.1033 

0.0482 

0.0000 

» 

4700 

O.OUOO 

0.0138 

0.0413 

0,089b 

0.151b 

0.213d 
1,2202 
2.34  !5 

0.2756 

0.4134  ) 

4600 

0.5512 

0.6890 

0 , 82o8 

0.9646 

1.1024 
2.204  7 

1.3780 

2. 4803 

1.5157 

4900 

1.6535 

1.7913 

1.9291 

2.0bi>9 

-0.1376 

5000 

0.0000 

-0.0344 

-0 . 0ob9 

-0, 1033 

-0.1240 

-0.13(9 

-0. 1 378 

5100 

-0.1419 

-0.1161 

-0.1516 

-0,151b 

-0.1516 

-0. 14ol 

-0.1419 

-0.1309 

5200 

-0.1033 

-0,075b 

-0,0482 

-0.02  7b 

-0.0138 

-0.0069 

0.0000 

5300 

1  1  7 

17 

5400 

-0,1000 

0,0000 

5500 

0,0000 

0.02 b2 

0.0524 

0,1047 

0,1571 

0.2094 

0.3142 

0.4189 

5b00 

0.5236 

0.6283 

0.6378 

1.0472 

1.2567 

1.3760 

1.7913 

2.2047 

5700 

2.7559 

0.2425 

5800 

0,0000 

0.0862 

0.1171 

0.1557 

0, 1H46 

0,2081 

0.2659 

5900 
6000 
b  1 00 

0.2797 

0.0000 

0.0000 

0.2880 

0.286b 

0,2673 

0.2343 

0,2136 

0,1406 

0,0799 

0.0262 

0.0524 

0.1047 

0.1571 

0.2094 

0.3142 

0,4189 

b200 
b  300 
6400 

0,523b 

2.7559 

0,0000 

o ,  e2b  3 

0.8378 

1,0472 

1 ,25e7 

1.3760 

.1.7913 

2.2047 

-0,04b2 

-  0 , 0  b  2  0 

-0,0758 

-0,0758 

-0,0690 

-0.0b20 

-0.0413 

6500 

-0,027b 

-0 , 02o7 

-O,0069 

0.00o9 

0,01 38 

0,0207 

0.027b 

0.0207 

b6O0 

O.OOOu 

" 

b700 

2 

6800 

1  1 

2  1 

-22,8700 

6900 

3  1 

2  2 

1  1 .5600  _ 

7000 

7100 

7200 

7300 

4  1 

1 

0.00 

1 

3  2 

15.1600^ 

, 

-  SLOT  DATA 

74uO 

0.125 

7500 

1.1483 

0.0833 

0,710 

7  bOO 

518.70 

1.20 

1.0 

7700 

0 

0 

0  . 

7600 

0 

0 

0 

7900 

8000 

8100 

0 

0 

0 

0 

0 

0 

0 

0 

0 

^ FLOW  PARAMETERS  j 

8  200 

0 

0 

0 

i 

8300 

0 

0 

0 

1 

8400 

0 

0 

0 

* 

6bO0 

1  l  ■ 

-  1  - 

- 

-  —  — 

— 

-  — 

■*■*>  —  -  - 

— 

6600 

THE  END 

ELEMENT  GEOMETRY 
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8600 
8700 
8800 
8900 
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9100 
9200 
9300 
9400 
9500 
9600 
9700 
9800 
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10000 
10100 
10200 
10300 
10400 
10500 
10600 
10/00 
10800 
10900 
11000 
11100 
11200 
11300 
11400 
11500 
1  looO 
11700 
11800 
119UU 
12000 
12100 
12200 

12300 
12400 
12500 
12600 
127  00 
12800 
12900 


7  00 
750 
800 
900 
1000 

1100 


1200 


1  oOO 
17  00 


: ^  ^iSlHARHjl'SuttJlKCK  bimS) 

■  rHJi  ( Ill  L2  L3  ) Sl  PAR  (  5  )  + 1 MULT+1  NCR  (LI )  *Z2 )  ^ 

lKfLi  EU  b)  PAR(L1,L2  L3  )  =  I PAR (6 )  +  (MULT*1NCR l LI  +NINCJ) 
IK(Ll*F0*7)  P  A  R  ( L 1  »  L  2  ,  L  3  ) =  1 P  A  K  (  7  )  +  (  M  U  LT  *  1 N  C  R  (  b  1 )  ■*  X  3  J 

IF  (  L3  #  £0  !  8  )  PAR IL 1  #  L2  r  b3 ) -IP AR  1 0 )  +  ( MIHiT*  J  NCR  (  b  1 )  ) 

TF  (IjI  V()  9)  PAR  (LI  »  L2  » L3 )-IPAR (9 )  +  (MULT* INCH  (Ll+NINC)  ) 
1FCL3 Ibo!  10)  PARUjl,b5,L3)alPAR(10)  +  (MULT*lNCR(bl  +  N£NU  ) 
It  (PAR(L1  ,1,2,10). LI. 0.0)  GOTO  500 
GOTO  750 

-CHECN^TU  PREVENT  DUPLICATION  OF  INITIAL  GEOMETRY 
IK  (L4 •  EU •  NPAR  )  GOTO  400 
PAk(Ll,L2,b3)  =  lPAi<U.i) 

CONTINUE 

CONTINUE 

CONTINUE 

WRJLTE^  (b  •  1100)  nUDNPtt,  1 1  .•  1  2  ,  N1NC,  WCASE  in- 

KORMAT  (///,7X,  'uEUMKTRY  MUUIFICaUON  #2X,12j//*/X,  II- 
I10,2X,'12=  '  ,110,2X, 'NINCr-  '  ,  12 , 2X  ,  NCALE-  »,12,//,9A, 

•  1  nCK  NUM  '  ,  3X  '  %  .STEP  • , 3X , '  DEGREE  '  ) 
l)U  1240  l=t,NlwC 

FORMAT  (1 2X , 1 2 , 7X , F5 . 2 , 4X , t  b .  3 ) 

CONTINUE  ia 

FORMAT1  (///?9X,  'PAKlJ  I  Ak^  'pAI<3  PAM  PA R5  P *\R6  PAR' 
'  PAR 8  PAR9  PAR  10',/,'  UPAR  '  ,10(11,  jX)) 

WRITE  (6,1300/  Ul  AR(i)r J=1,10J 
FORMAT  ('  1N.GE.  '  ,  10(.Fb.2)  ) 

WHILE  (6,1400)  , 

FORMAT  ('  IN  C8  ') 

DO  1700  1,1  =  1  ,  NJ  DC 
DU  1600  E2= 1 , uCanE 

L3= ( L 1  - 1  )  ♦  NCAJSE+L2  ( 

WRITE  (6,1500)  L1,L3,(PAR(L1,L2,J),0=1»10)  j 

FORMAT  (/,X,I2,X,l2,X,10(to.2  ,/,9X,  EL-  ,oA, 

•CU=  *,6X,'CL/CD=  • ,6X,'Cn-  ) 

CONTINUE 

S TOP MODI F 1  CATION  COMPLETE'  I 

END 


AFPKgDIX  ( iii)  RAl'TO^I  ■  "TITTO  C-1?1'F.*AW 


D.4.9  RANDU  Subroutine 

The  RANDU  subroutine  computes  a  pseudo-random  number,  as  a  single 
precision  value  uniformly  distributed  in  the  range: 

0.0  .LB.  value  .LT.  1.0 

Format: 

CALL  RANDU (il,i2,x) 

Arguments: 

il,i2 

INTEGER*2  variables  or  array  elements  that  contain  the  seed  for 
computing  the  random  number. 

x 

a  real  variable  or  array  element  where  the  computed  random  number 
is  stored. 

Notes: 

1.  The  values  of  il  and  i2  are  updated  during  the  computation  to 
contain  the  updated  seed. 

2.  The  algorithm  for  computing  the  random  number  value  is  as 
follows: 

If  21=0 ,  sot  generator  base 

X (n+1)  =  2**16  +  3 
otherwise 

X (n+1)  =  (2** 16+3)  *  X (n) mod  2**32 
Store  generator  base  X(n+1)  in  11,12. 

Result  is  X ( n+1 )  scaled  to  a  real  value  Y(n+1),  for  0.0  .LE. 
Y (n+1)  .LT.  1. 
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